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ABSTRACT 

We report on optical and X-ray observations of the accretion powered ms pulsar IGR J00291+5934 
in quiescence. Time resolved /-band photometry has been obtained with the 4.2 m William Herschel 
Telescope, while a 3 ks Chandra observation provided contemporaneous X-ray coverage. We found 
an unabsorbed 0.5-10 keV X-ray flux of lxl0~ 13 erg cm~ 2 s _1 which implies that the source was 
in quiescence at the time of the optical observations. Nevertheless, the optical /-band light curve of 
IGR J00291+5934 shows evidence for strong flaring. After removal of the strongest flares, we find 
evidence for an orbital modulation in the phase folded /-band light curve. The overall modulation 
can be described by effects resulting from the presence of a superhump. Comparing our lightcurve 
with that reported recently we find evidence for a change in the quiescent base level. Similar changes 
have now been reported for 4 soft X-ray transients implying that they may be a common feature of 
such systems in quiescence. Furthermore, the maximum in our folded lightcurve occurs at a different 
phase than observed before. 

Subject headings: stars: individual (IGR J00291+5934) — accretion: accretion discs — binaries: close 
— stars: neutron — X-rays: binaries 



1. INTRODUCTION 

Low-mass X-ray binaries (LMXBs) consist of a com- 
pact object, either a neutron star or a black hole, that 
accretes from a late type companion star. The compan- 
ion star typically has a mass < 1 M . Stellar evolu- 
tion theory predicts that neutron star LMXBs are the 
predecessors of the recycled millisecond radio pulsars 
(jRadhakrishnan fe Srinivasanl [1981 lAlpar et all Il982f ). 
This link between the neutron star LMXBs, their tran- 
sient cousins, the X-ray transients, and the millisecond 
radio pulsars has been established by the discovery of the 
first accretion powered millisecond X-ray pulsar in 1998 
(IWiinands fc van der Klisl [l998; Cha krabartv fc Morgan! 
1998). 

Subsequently, more of these systems have been found. 
At the time of writing there are ten accretion pow- 
ered millisecond X-ray pulsars known. In three of 
these sources pulsations are detected only intermit- 
tently (the transients HETE J1900. 1-2455, Aql X- 
1 and SA X J1748.9-2021 in the globular cluster 
NGC 6440; iGallowav et all [20071 iCasella et all [20071 . 
lAltamirano et al.ll2007[) . To date all o f them are found in 
X-ray transients ( see IWiinan ds 2006 for an observational 
overview). These transients flare-up during episodes of 
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enhanced mass accretion rate onto the neutron star. It is 
thought that during outburst the enhanced mass accre- 
tion rate suppresses the radio pulsar mechanism preclud- 
ing the detection of pulsed radio emission. Radio emis- 
sion of these transients has been detected in outburst, 
but this is associated with synchro tron emission from a 
radio jet (e.g. S AX J1808.4-3658. iRupen et all l2002ri 
XTE J0929- 314lRupen et al.ll2002at IGR J00921+5934. 
lPoolevll2004l . iFender et al.ll2004h . When the mass accre- 
tion rate stops/becomes very low the radio pulsar mecha- 
nism should turn on. However, despite deep searches us- 
ing sensitive radio telescopes at epochs when one of these 
transients, SAX J1808. 4-3658, was in qui escence no ra- 
dio m illisecond pulsations have been found (jBurgav et al.l 
2003). The non-detection has been ascribed to the pres- 
ence of absorbing material close to the binary. 

In quiescence these systems are very faint in the 
optical. The counterpart can often not be detected 
(e.g. XTE J0929-314 and XTE J1814-3 31, the latter 
sourc e was not detected down to /?=23.3; iKrauss et alJ 
I2005f ). For the system XTE J 175 1-305 the optical coun- 
terpart wBS_not_chscoyered in outburst nor in quies- 
cence (jJonker et al.l l2003| ) . T here are two notic eable 
exceptions; SAX J1 808.4-3658 dHomer et al.l 1200 lft and 
IGR J00291+5934 (TFbx fc Kulkarnil 120041: iRoelofs et al.l 
120041 : ISteeghs et al.l l2004t iBikmaev et all 120051 ). For a 
comprehensive overview of the outburst and initial quies- 
cence o bservations of IGR J002 91+5934 see lTorres et al.l 
(|2007f) . iD'Avanzo et all (|2007[ ) recently reported multi- 
band quiescent optical and near-infrared observations of 
IGR J00291+5934. 
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Using optical observations of SAX J1808. 4-3658 in qui- 
escence, Homer et al. (2001) found evidence for a 9-15 
per cent semi-amplitude modulation (the observed am- 
plitude depends on the photometric band that is used 
for the observations). Those authors proposed that this 
would be due to X-ray irra diation from the neut ron star. 
However, as pointed out bv lBurderi et al.l (|2003[ ), assum- 
ing that SAX J1808. 4-3658 was indeed in quiescence at 
the time of the optical observations, the quiescent X-ray 
flux of SAX J1808.4-3658 is too low by two orders of 
magnitude to explain the modulation in terms of X-ray 
heating. The absence of a double-humped morphology 
rules out that the modulat ion is due to ellipsoidal varia- 
tions. iBurderi et al.l (|2003f ) proposed that the irradiation 
is caused by a turned-on radio pulsar instead of the qui- 
escent X-ray emission. Evidence for a turned-on radio 
pulsa r, albeit indirect in this case (see ICampana et al.1 
l2004f) . would reinforce the link between the LMXBs and 
the millisecond radio pulsars. Furthermore, a turned-on 
radio pulsar would have an important effect on the evo- 
lution of the mass-losing d onor star, altering the evolu- 
tionary path of the binary (| Ruder man et al.l Tl989). Due 
to the absence of pointed X-ra y observations at the time 
of the optical observa tions of iHomer et alj (|2001h and 
ICampana et al.l (|2004f ). low level X-ray activity could 
have remained unnoticed in X-rays but would heat the 
side of the companion star facing the neutron star, pro- 
viding an alternative explanation for the observed op- 
tical modulations. However, simultaneous X-ray and 
optical observation s of SAX J1808. 4-3658 reported by 
Hci nke et al.l (|2007f ) forego this possibility. 

In this Manuscript, we present phase resolved pho- 
tometric observations of the optical counterpart of 
ICR J00291+5934 in quiescence obtained with the 4.2 m 
William Herschel Telescope (WHT). In addition we 
present our analysis of a short contemporaneous Chan- 
dra observation of IGR J00291+5934 obtained with the 
aim to determine the contemporaneous X-ray flux. 

2. OBSERVATIONS, ANALYSIS, AND RESULTS 

2.1. Chandra X-ray observations 

We have observed IGR J00291+5934 with the back- 
illuminated S3 CCD-chip of the Advanced CCD Imag- 
ing Spectrometer (ACIS) detector on board the Chandra 
satellite. The observation started on Sept. 13, 2006 at 
13:33 (UTC; MJD 53991). The data telemetry mode was 
set to very faint to allow for a thorough background sub- 
traction. Due to windowing of the ACIS-S CCD a frame 
time of 0.4104 s has been used. We have reprocessed and 
analysed the data using the CI AO 3.4 software devel- 
oped by the Chandra X-ray Center using CALDB ver- 
sion 3.3.0.1 to take full advantage of the very faint data 
mode. In our analysis we have selected events only if 
their energy falls in the 0.3-7 keV range. The 0.3-7 keV 
background count rate was always lower than 0.4 counts 
s _1 during our observation. The net on-source exposure 
time is 2.88 ks. 

Using wavdetect we detected two sources 
CXC J002911+593420 and IGR J00291+5934. We 
detect 22 photons from IGR J00291+5934 during 
the 2.88 ks observation, leading to a source count 
rate of 7.6xl0~ 3 cou nts s _1 . Using xspec version 
11.3.2p {A rnaudl [i~996t ) we have fitted t he spectrum of 
IGR J00291+5934 using Cash statistics (|Cashlll979ft to 




Fig. 1.— Optical /-band finder chart of IGR J00291+5934. 
North is up East to the left, the scale of the image is indicated 
in the top-right part of the figure. The optical counterpart to 
IGR J00291+5934 is indicated with two vertical lines. The stars 
labeled CI to C4 are comparison stars with /-band magnitudes 
18.41±0.02, 18.55±0.03, 19.34±0.05 and 19.41±0.03, respectively. 
The star marked with C5 is a nearby comparison star of similar 
brightness as IGR J00291+5934 (I = 21.80 ± 0.05). 

an absorbed powerlaw model. We held fixed the inter- 
stellar extinction to the value of 4.6xl0 21 cm~ 2 favored 
by Torres et al. (2007) and the powerlaw index to 2. We 
derive an absorbed 0.5-10 keV source flux of 7xl0~ 14 
erg cm" 2 s _1 and an unabsorbed 0.5-10 keV flux of 
lxlO -13 erg cm~ 2 s _1 . The source flux is consi stent 
with that derived previously by iJonker et alj (|2005l ) and 
Torres et al. (2007). We searched for variability in the 
rate of arrival of the photons but we fo und none. A 
Kolmogorov-Smirnov (jPress et al.l Il992h test showed 
that the probability that the data are consistent with 
the null-hypothesis of a constant photon arrival rate is 
63%. We conclude that the source was in quiescence 
during our optical observations. 

2.2. WHT optical observations 

We obtained Harris /-band images using the Auxiliary 
Port Imager (AUX) instrument mounted on the 4.2 m 
WHT telescope at the Roque de Los Muchachos Obser- 
vatory, La Palma, Spain. On September 13 and 14, 2006 
(MJD 53991 and 53992 UTC) observations with an ex- 
posure time of 600 s, totalling 10.7 hours of data, were 
obtained. The observing conditions were good with a 
photometric sky and a seeing varying between 0'/7-l'.'6 
and between O'.'7-l'.'O and with a mean seeing of l'.'O and 
/ /74 on Sept. 13 and Sept. 14, respectively. Given the 
faintness of the quiescent counterpart we have rejected 
images with seeing >1". This left us with 24 images for 
Sept. 13 and 40 for Sept. 14. We used a 2 x 2 on-chip 
binning, therefore AUX delivered a circular field of view 
with a diameter of 1.'8 sampled at O'.'22/pixel. Frames 
were debiased and then flatfielded using dome flatfield 
observations. 

We applied both aperture photometry and point- 
spread function fitting on each of the 600 s images us- 
ing DAOPHOT in IRAF 3 to compute the instrumental 
magnitudes of the detected stars. The results obtained 
with both methods are consistent with each other. Flux 

3 IRAF is distributed by the National Optical Astronomy Obser- 
vatories 
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Phase folded lightcurve 

Fig. 2.— Optical /-band folded light curve of IGR J00291+5934 
(stars) obtained with the 4.2 m WHT on Sept. 13 and 14, 2006. 
Strong flaring of IGR J00291+5934 can be seen. The strong flaring 
precludes the detection of clear trends in the /-band magnitude as 
a function of orbital phase. The open circles shows the /-band 
magnitude of a nearby star of comparable brightness (star C5 in 
Figure [TJ measured simultaneously. The points for star C5 have 
been scaled down by 0.5 magnitude for display purposes. Phase 
zero corresponds to the epoch of 90° mean longitude as given by 
Galloway et al. (2005). 

calibration of the field was performed by observing two 
Landolt standard stars and differential photometry was 
used to derive the source flux variability as a function 
of time. The photometric results given here are with re- 
spect to the four field stars shown in Figure [T] These are 
the brightest isolated stars available in the unvignetted 
field of view of AUX that were recorded in the linear 
regime of the CCD. 

We plot the /-band magnitude of IGR J00291+5934 
(crosses) and those of a nearby comparison star of similar 
brightness (open circles; star C5) as observed on Sept. 13 
and 14, 2006 in Figure [2] folded on the orbital period of 
I GR J00291+5934 We have used the pulsar ephemeris 
of iGallowav etaD (f2005h . The error in propagating the 
orbital ephemeris to the time of our observations is «0.01 
in phase. Phase zero is superior conjunction of the neu- 
tron star (i.e. the epoch of 90° mean longitude). The 
brightness of the comparison star (star C5 on Figure [1} 
is consistent with being constant. In contrast, the folded 
lightcurve of IGR J00291+5934 displays several large 
flares (see also Figure |3|) . The average /-band magni- 
tude of IGR J00291+5934 and the variance therein is 
21.83 and 0.18 magnitudes, respectively. 

In an attempt to determine whether there are varia- 
tions in the quiescent optical brightness due to a change 
in aspect of the Roche-lobe filling companion star as a 
function of the orbital phase, we have removed flares from 
the lightcurve of both nights. In order to come to a defi- 
nition of a flare we have varied the magnitude threshold 
above which we identify a data point as a flare. Next, for 
a range of potential magnitude thresholds we fold and av- 
erage the data in 10 orbital phase bins and we determine 
the x 2 and the number of degrees of freedom (d.o.f.) of 
a fit of a sinusoid plus a constant to the folded data. We 
show the resuls of this in Figured] It is clear that the x 2 
strongly decreases when more and more stringent mag- 
nitude thresholds are taken. The reduction in \ 2 levels 
off around a threshold magnitude of 21.75 (see Table [1]). 



Fig. 3. — The image corresponding to the strongest flare visible 
in Figure [2]is displayed on the left. On the right is the image taken 
20 minutes after the flare. North is up East to the left. 



TABLE 1 
The effect of the definition of a 
flare on a fit of a constant plus 
sinusoid to the folded lightcurve. 



Magnitude limit 




X 2 


D.o.f. 


21.95 


1.1 


13.2 


12 


21.90 


1.77 


35.4 


20 


21.85 


2.44 


73.2 


30 


21.825 


2.82 


98.7 


35 


21.80 


3.05 


134 


44 


21.775 


3.23 


152 


47 


21.75 


3.29 


161 


49 


21.725 


3.43 


175 


51 


21.70 


3.56 


185 


52 


21.65 


4.23 


224 


53 


21.60 


4.95 


272 


55 


21.55 


5.82 


326 


56 



For more stringent magnitude limits the decrease in \ 2 
is proportional to the reduction in the number of degrees 
of freedom as indicated by the straight dashed line in 
Figure |4] Such behavior is expected when clipping data 
points that are not flares. 

Using 21.75 as our magnitude threshold above which 
we define data points as flares, the average J-band mag- 
nitude and the variance therein becomes 21.90 and 0.09 
magnitudes, respectively. The variance is slightly higher 
than ~0.05 which is found for several stars of similar 
brightness that do not vary, implying that variability due 
to intrinsic variations is still present in the lightcurve of 
IGR J00291+5934. The resultant folded light curve is 
shown in Figure [5J The most striking feature is the pres- 
ence of a clear orbital modulation. In order to quantify 
the modulation we have fit the folded lightcurve with a 
fit function consisting of a constant plus a sinusoid. An 
F-test gives that the improvement of a fit of a constant 
plus sinusoid with respect to a constant has a probability 
of 2% to occur due to chance. We derive (6 ± 1) x 10~ 2 
magnitudes and 0.34± 0.03 for the semi-amplitude and 
phase of the sinusoid, respectively. The value of the con- 
stant is 21.91±0.01 magnitudes. In addition to the si- 
nusoidal variation there is suggestive evidence for an in- 
crease in source brightness at phase 0.9. In principle, an 
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Fig. 4. — The x 2 of the fit of a constant plus sinusoid to the phase 
folded light curve of IGR J00291+5934 plotted against the number 
of degrees of freedom (d.o.f.) in the fit for the various magnitude 
thresholds that we investigated (see Table [1} . The drawn line is 
the linear correlation described by the data points below d.o.f.=50. 
The transition between the linear correlation and a more steep 
correlation occurs around magnitude 21.75. 
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Fig. 5. — Phase folded optical /-band light curve of 
IGR J00291+5934 after excluding flares in the light curve (i.e. data 
points with / <21.75). The dashed line is the best— fitting model 
of a sinusoid plus a constant. The sinusoid reaches a maximum at 
phase 0.34 ± 0.03. For clarity two orbital cycles are plotted. 

alignment of small flares in the phase bin centered on 
phase 0.9 can have caused an increase above the sinu- 
soidal modulation. Ellipsoidal modulations due to the 
Roche-lobe filling companion star are not detected. 

3. DISCUSSION 

Comparing the source X-ray flux during a 3 ks-long 
Chandra observation of IGR J0029f +5934 with that ob- 
served on previous occasions when the source was in 
quiescence, we conclude that IGR J00291+5934 was 



in quiescence on Sept. 13, 2006. Optical /-band ob- 
servations acquired with the 4.2 m WHT telescope on 
Sept. 13 and 14, 2006 reveal strong variability with 
flares of ~1 magnitude on timescales of tens of min- 
utes. After removal of flares we have phase folded our 
data. A clear sinusoidal variation is present. In our 
data set the orbital phase at w hich the source bright - 
ness is maximal is 0.34 ± 0.03. iD'Avanzo et all (|2007l ) 
report a maximum consistent with phase 0.5. Further- 
more, their average /-band magnitude of 22.4±0.2 for 
IGR J0029I+5934 is 0.50±0.22 lower than our average 
magnitude (after removal of flares larger than 7=21.75). 
Hence, th e brightening is only s ignificant at the 2.3 a 
level and ID'Avanzo et alj (|2007f ) obtained the /-band 
observat ions under non- p hotom etric conditions. Inter- 
estingly, ID'Avanzo et al.l (|2007f ) do not report evidence 
for flaring. The lack of flaring together with the lower 
average brightness suggests that the source showed ac- 
tivity at the time of our optical observations above that 
found bv ID'Avanzo et all <|2007h . 

In order to investigate the nature of the flares and 
the observed variability in the phase folded optical 
light curve, we compared the quiescent properties of 
IGR J00291+5934 observed by us with tho s e rep orted 
by Torres et al. (2007) and ID'Avanzo et~aT] (|2007h and 
with those report ed for several other sources, most no- 



tably Cen X-4 (jChevalier et alJ I1989L iShahbaz et "all 
1993). Optical flares similar to those observed in 
IGR J00291+5934 in quiescence have been reported for 
sever a l sources (e.g. | Chevalier et al.l I1989L IZurita et al.l 
120031 iHvnes et al.1 120031 )! However, flares as large as 
the largest that we observe in IGR J00291+5934 have 
not been reported to our knowledge. In addition, we 
note that a variation of about half a magnitude in 
the average brightness similar to that mentioned above 
for IGR J00291+5934 has been reported for Cen X- 
4, XTE J21 23-058 and the black hole X-ray transient 
A 0620-00 (IChevalier et all I1989L iTomsick et all 12004 
ICantrell et al.ll2008l . respectively) . 

Strong flares of similar duration (tens of minutes) 
have been observed in m id and late M dwarfs (e.g. see 
Rockcnfcl ler et all 120061 and references therein) . How- 
ever. IHvnes et all ( 20031 ) showed that the flares in qui- 
escent LMXBs are associated with the accretion disk. 
Indeed, the combination of strong flaring and the lack of 
clear ellipsiodal varia tions in the phase folde d light curve 
(this work as well as ID'Avanzo et aT1l2007f ). imply that 
the quiescent light is not dominated by a Roche lobe fill- 
ing, but otherwise unperturbed donor star. If the flares 
are from the donor star their amplitude must be larger 
still. Instead, (superhumps from) a residual accretion 
disk, the accretion stream, and/or effects from an irradi- 
ated donor star are potentially important contributors to 
the quiescent /-band light. Torres et al. (2007) derive a 
similar conclusion on the basis of the quiescent intrinsic 
R — K color and constraints on the donor star from the 
pulsar mass function. 

The absence of ellipsoidal variations in the phase 
folded light curves of both IGR J00291+5934 as well as 
SAX J1808. 4-3658 is consistent with the very low contri- 
bution of a (non-irradiated) brown dwarf to the optical 
light. Overall, the phase folded quiescent /-band light 
curve of IGR J00291+5934 that we find resembles that 
of the average quiescent V^-band light curve of Cen X- 
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1989) in 1984-1988 (see 



1989) show that a model 



4 observed by I Chevalier et al.l 
their figure 7). IChevalier et al.l 
where a large fraction of the optical light comes from 
the heated hemi-sphere of the companion star or an ac- 
cretion disk attenuated by electron scattering in an ac- 
cretion wake can describe the data well. In particular 
such a model can explain the lower brightness at orbital 
phases 0.5-0.7. The heating of the companion star in 
ms X-ray pulsars in quiescence has been ascribed to th e 
turn-on of an active radio pulsar (|Burderi et al.ll2003f) . 
Such an effect might be underlaying the folded lightcurve 
of IGR J00291+5934, however, due to the (likely) higher 
level of activity during our observations additional ef- 
fects seem to be important. An important clue comes 
from the phase difference between the maximum of the si - 
nusoidal modulation reported by iD'Avanzo et al.l (|2007f ) 
and that reported in this work. This suggests the pres- 
ence of a modulation that has a period unequal to the 
orbital period, such as a superhump. Potentially, the 
overall morphology of our folded lightcurve can indeed 
be described by effects originating in a superhump. Su- 
perhumps have been detected in the black h ole candidate 
XTE Jl 118+480 while near quiescence (IZurita et ail 
120021) . As explained in IWhitehurst k Kind (|1991h. su- 
perhumps can occur in systems when the ratio between 
the companion star and the compact object mass is < 
0.3. An accretion stream that perturbs the outer disk 
facilitates the excitation of ti dal torques that are settin g 
off the superhump resonance ( Whitch ursTfc KindH99lh . 



The hint for a brightening of the source around phase 0.9 
visible in Figure [5] can be explained by the presence of 
an accretion stream and hot-spot while the mass ratio is 
very likely (much) lower than 0.3. However, from studies 
of Cataclysmic Variables it is known that the superhump 
only develops during outburst, most likely since the mass 
and hence angular momentum transfer rate through the 
disk in quiescence is too low to allow for a large extent of 
the disk. Hence, it is not clear whether all the conditions 
necessary for the occurrence of superhumps are fullfilled 
in IG R J00291+5934 durin g our observations. Interest- 
ingly, iNeilsen et alJ (|2007f ) also found evidence for the 
presence of an eccentric, precessing disk in the black hole 
candidate A 0620-000 in quiescence. As superhumps can 
potentially explain the overall variability observed in the 
optical in several systems while in quiescence, further in- 
vestigation is necessary to test whether during phases of 
(enhanced) mass transfer in quiescence a superhump can 
develop in quiescent LMXBs with mass ratios less than 
0.3. 
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